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The oscillation of tunnel splitting in the biaxial spin systems with the magnetic field 
along the hard anisotropy axis is analyzed within the particle mapping apporach, 
rather than in the (9, 4>) spin coherent-state representation. In our mapping pro- 
cedure, the spin system is transformed into a particle moving in the restricted S 1 
geometry whose wave function subjects to the boundary condition involving addi- 
tional phase shift. We obtain the new topological phase that plays the same role as 
the Wess-Zumino action in spin coherent-state representation. Considering the inter- 
ference of two possible trajectories, instanton and anti-instanton, we get the identical 
condition for the field at which tunneling is quenched, with the previous result within 
spin coherent-state representation. 
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For the last decade, there have been a lot of works on the macroscopic magnetization 
reversal in nanomagnetic single-domain particles with many spins both theoretically!'! and 
experimentallyiHi. It is well-known that, at low temperature, the rate of such reversal 
is dominated by quantum tunneling through the barrier, while at high temperatures, it is 
governed by classical thermal activation according to the law exp(—E /kBT), where E is the 
barrier energy and T is the temperature. Moreover, two kinds of tunneling mechanism exist 
in tunneling regime, the macroscopic quantum tunneling (MQT) from a metastable state to 
a stable one, and the macroscopic quantum coherence (MQC) between the two equivalent 
ground states. These have been treated by spin coherent-state path integral formalismi or 
direct mapping onto a particleEHl in uniaxial and biaxial spin systems. 

Here, we focus on the MQC of biaxial spin system like Fes which shows the exotic 
phenomena. In the absence of magnetic field, Loss et. al.0 and Delft et. al.0 have shown 
that the tunneling is quenched for half-integer spins, and allowed for integer spins using 
the coherent state path integral. Such spin parity effect is orginated from the topological 
Wess-Zumino phase in Euclidean action in path integrally! and is directly related to the well- 
known Kramer's degeneracy. Later, Garg0 first have shown, also within the coherent-state 
path integral, that in the presence of magnetic field along the hard anisotropy axis, i.e., in 
the system decribed by Hamiltonian 

n = hs 2 z + k 2 s 2 x - hs z , (i) 

where k\ and k 2 are positive anisotropy constant with k\ > k 2 , the tunneling spilitting also 
vanishes when the external field H takes the values 



H = 2k 2 Vl- A(5 -n-1/2), (2) 

where n = 0, 1, ■ • ■ , 25 — 1, A = k 2 /k\ and 5 is the total spin. Very recently, this system 
has been conformed by experimental, and also analyzed in more detailS. These oscillations 
of level splittings are also due to the topological phase factor in spin coherent-state path 
integrals. These seemingly simple level crossing can be interpreted as the interference of two 
different classical trajectories in Euclidean space. 



Another approach, the direct mapping onto a particle, is also useful in that this enables us 
to treat the spin system more quantum mechanically, as well as gives a constant mass in the 
biaxial spin system. In this paper, we derive Eq. (^j using the direct particle mapping rather 
than coherent-state path integral. In this procedure, the boundary condition of particle wave 
function in the presence of external field along the hard axis is modified by a topological 
phase shift in comparison to the zero-field case. We can interpret that a particle moves on a 
restricted S l geometry with wave function having not the periodic boundary condition but 
the additional geometrical phase. From this boundary condition, we can derive the overlap 
integral which produces the classical action in Euclidean space. Considering the interference 
of two possible trajactories, instanton and anti-instanton, we find the field values at which 
the tunneling vanishes. These values are just the same as the previous results obtained 
within the spin coherent-state path integral formalism. 

Consider first the simple uniaxial spin system which has hard axis anisotropy with mag- 
netic field applied along the hard axis, in order to get an idea in a simpler system. The 
corresponding spin Hamiltonian is expressed as: 



Note that this system does not correspond to tunneling. This Hamiltonian has eigenvalues 
e m = m 2 — Hm in ^-representation, where m = —S, —S — 1, ■ ■ ■ , S — 1, S, and S is the total 
spin number. The eigenvalues have many level crossings as the field increases. Although we 
already knew all the energy levels, for the pedagogical purpose, we perform the mapping 
onto a particle by defining the characteristic function 



Following the well-known mapping procedure, we obtain the Schrodinger-like equation 



T~C = S 2 — HS Z . 



(3) 




(4) 



d 2 V H 2 



(5) 



d(j) 2 4 



where 




(6) 
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The wave function \1/ satisfies the boundary condition for 2-7r-rotation 

+ 2tt) = ±V{<f>) exp(i7iH) (7) 

where plus sign corresponds to bosons, while minus sign to fermions. This equation corre- 
sponds to the motion of particle with mass m = 1/2 moving in the potential V(4>) = H 2 /4 
on the constrained geometry S 1 . In above example, the phase factor appearing in the wave 
function, Eq. ([]), plays a crucial role in yielding the original eigenvalues, e m = m? — Hm. 
Although the level crossing is merely the quantum mechanical property in spin systems, 
such effect can be interpreted as the result of the geometrical phase in the wave function of 
particle. 

Another uniaxial spin system, Mn^Ac, whose Hamiltonian is given by 7i = —S 2 — HS Z 
can be considered. Following the above procedure, we obtain the equation describing the 
particle with negative mass. This corresponds to the motion in the inverted potential having 
(2S + l)-eigenvalues — e m but in opposite order to those of original spin system. The level 
crossing in this case can also be shown by the same method. These observations are to be 
extended to the biaxial spin systems with external magnetic field applied along the hard 
anisotropic axis. 

Let's consider the biaxial spin Hamiltonian with the magnetic field in the hard ^-direction 

H = -DSl + ES 2 Z - HS Z} (8) 

where D and E are easy- and hard-axis anisotropy constants satisfying D, E > 0. This 
system also can be mapped onto a particle system using the usual mapping procedurei. We 
define the characteristic function 

S imcf> 

$(0)= E (9) 
m =s y(S- m)\(S + m)\ 

where c m is the wave function of spin system. Then, we obtain the differential equation for 
$: 

- P (<f>)-7-2 + QW-il + = ^ ( 10 ) 



where 

P(0) = 1 - A + Asin 2 0, 

Q((f)) = 2iSh + \(s- -)sin20, (11) 
= -^(S 2 + S) - ^(S 2 - S) cos 20, 

with A = D/(D+E), h = H/2(D+E)S, and e is the reduced energy eigenvalue. To eliminate 
the first-order derivative of $ and make the particle mass constant, we can introduce the 
new wave function 

W) = $(0)exp(-G(0)), (12) 

and change the variable from to x = x(<f)), so that the Eq. (pi]) can be transformed as 

d 2 m 



dx 2 



+ V(x)ty = e-V, (13) 



where 



d6 

1=1 /i x I x • TZ = sn-^sin^V^), (14) 
VI - A + Asur # 

Here, sn(x, k) is the usual Jacobian Elliptic funtion with period of 4/C(A;), where /C is the 
complete elliptic function of first kind. Due to Eq. (|14"D, sin0 and cos0 can be changed into 
snx and cnx, respectively. The resulting effective potential and wave function are given by 

V(x) = 1 - „ | X 2 S(S + l)sn 2 xcn 2 a; - h 2 S 2 

1 — A + Asmx \ 

+2i\hS (^S + snxcnx^ — A(S' 2 + S')cn 2 x 



\ft(x) = $(x)(l - A + Asn 2 x) 2 

dnx'cir' 



x exp ^ — J — 



(15) 



(16) 



A + Xsn 2 x' J 

Note that the effective particle potential corresponds to the complex potential. The real 
part of V(x) is the periodic potential with period 4/C. The imaginary part of V(x) is related 
to the topological phase shift in wave function in the manner of Aharonov-Bohm effect, but 
does not contribute to the tunneling rate, which will be clear later. Since the integral 



X + 4!C dnx'dx' 



2tt 



1-A + AsnV v^T - T 
wave function \1/ satiesfies the following boundary condition: 

h 



(17) 



W(x + 4JC) = ^(x)exp 



2mS 1 



(18) 



Now, we consider the path integration in Euclidean space in order to get the tunneling 
rate. The amplitude from the state | x > at Euclidean time r = Tq to the state | x' > at r is 



< x',r \ x, r >= / V[x(t)] Y[ < Xk | > exp ( — / drTi 

J k=1 ' V Jr 



(19) 



where x® = x and xn = The overlap integral between kth and (k — l)th space elements 
can be written in small x^ — Xk-\ limit as 

2nS , , / /i ~ 



< £fc | ^fc-i >= exp 



a; fc - x fc _i) 1 



(20) 



Therefore, the product of overlap integral between fcth and (k — l)th space elements in Eq. 



(|19|) is given by 



N 



Y[ < x k \ x fc _i > = exp 
fc=i 

= exp 



-If? 6 * I 1 



h 



h 



-i—r^ ( drx f 1 . 

2/C K \ ^T^rj 



(21) 



where e = — Tfc_i. The imaginary-time effective action which gives the tunneling rate is 
therefore written as 



7T 



S 



S R = -i— 1 



h 



J dx + J dr[4x 2 + V(x)}. 



(22) 



2/C VT^x, 

We consider the interferences of two possible trajectories of instanton and anti-instanton 
which move from x = to x = 2/C and to —2/C in the real part of inverted potential. In this 
case, the real parts of tunneling rate are factored out and we only focus on the imaginary 
part of effective action in Eq. 



ImSV 



nS 



1 



h 



J dx + J drlmV(x). 



(23) 



2/C V VT^X, 

Here, the contribution of second integral, that is, imaginary potential, is zero, due to its 
oddness. Therefore, we obtain the interferences of two trajectories which vanishes when 
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h = VT^\ S n 1/2 , n = 0, 1, • • • , 2S - 1. (24) 

which is exactly same as those in the previous calculations within (8, <fi) spin coherent-states 
represent ationBlll 

Here, we note that the variable 0, and therefore x seems be an azimuthal angle. When we 
choose the Hamiltonian as Eq. (||) among the possible Hamiltonians that have same physics, 
can be interpreted as the usual azimuthal angle in the ^-representation, since z-axis is 
a hard axis. But in this case, while the particle mass is positive, the effective potential 
becomes complex, and therefore, the energy eigenvalues are also complex. However, the 



boundary condition with the additional phase, Eq. (|18|) just exactly singles out the real 
energy eigenvaluesi. We can choose the alternative Hamiltonian other than Eq. (§) in order 
for the imaginary potential not to appear in Schrodinger equation. For example, let us 
consider the Hamiltonian 

n = -AS 2 - BS 2 X - HS y , (25) 

where A and B are positive constants with A > B > 0. We can also make a direct mapping 
onto a particle as above. In this case, the variable <fi in the characteristic function does 
not play a role as an azimuthal angle, so that the tunneling picture in original spin system 
does not have one-to-one correspondence with an effective particle systenfi However, the 
effective eigenvalue equation is given by 

d 2 V(x) 



dx 2 



- V{x)V{x) = e*(x), (26) 



where 



,2. 



cn"x 

V(x) 



sn 2 x 



b 2 S(S + l)sn 2 x - 2bhS ( S + l - ) snx + h 2 S 2 



+bS(S + l)cn 2 x + hS(2S + l)snx, (27) 
with b = B/A and h = 2AS. The boundary condition of particel wave function is given 

V(x + 4JC) = ±*(x). (28) 
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As shown above, the Hamiltonian gives the real potential with negative mass. Furthermore, 
the field- dependent additional geometrical phase in the boundary condition of wave function 
does not appear by this Hamiltonian. As discussed in the uniaxial case, this can be regarded 
as the motion in the inverted potential with eigenvalues having opposite sign and in opposite 
order to the original spin system. Through the simple numerical calculation, we can obtain 
the (25 + 1) inverted energy eigenvalues, which show the oscillation of tunnel splitting with 
external field. This is represented in Figure. [I] 

In conclusion, in this paper we consider the oscillation of tunnel splitting in spin system 
within direct particle mapping approach. In simple uniaxial model, the eigenvalues showing 
the level crossing with field can also be obtained from the topological phase in particle wave 
function occuring when a particle completely moves around a unit circle. For biaxial spin 
system with magnetic field along the hard axis, particle wave function has the topological 
phase term due to field in the boundary condition as in the uniaxial case, and this phase 
gives the imaginary part of Euclidean action analogous to the Wess-Zumino phase in (6, <ft) 
spin coherent-state representation. Therefore, in particle mapping formalism the oscillation 
of tunnel splitting can be understood as the topological effect of a particle moving in a 
periodic potential, that is, the constrained S 1 geometry with the boundary condition giving 
the geometrical phase. 
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FIG. 1. The lowest (2S + 1) energy eigenvalues of Eq. (p6|). This structure is exactly same 
as that of original spin system expressed as Hamiltonian, Eq. ( p5[ ) in inverted order. The points 
where the solid and dotted lines meet represent the level crossing. 
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Fig. 1 



